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We  demonstrate  the  use  of  an  optical  frequency  comb  to  coherently  control  and  entangle  atomic  qubits. 

A  train  of  off-resonant  ultrafast  laser  pulses  is  used  to  efficiently  and  coherently  transfer  population 
between  electronic  and  vibrational  states  of  trapped  atomic  ions  and  implement  an  entangling  quantum 
logic  gate  with  high  hdelity.  This  technique  can  be  extended  to  the  high  held  regime  where  operations  can 
be  performed  faster  than  the  trap  frequency.  This  general  approach  can  be  applied  to  more  complex 
quantum  systems,  such  as  large  collections  of  interacting  atoms  or  molecules. 
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The  optical  frequency  comb  generated  from  an  ultrafast 
laser  pulse  train  has  revolutionized  optical  frequency  met¬ 
rology  [1-4]  and  is  now  playing  an  important  role  in  high 
resolution  spectroscopy  [5].  The  spectral  purity  yet  large 
bandwidth  of  optical  frequency  combs  also  provides  a 
means  for  the  precise  control  of  generic  quantum  systems, 
with  examples  such  as  the  quantum  control  of  multilevel 
atomic  systems  [6,7],  laser  cooling  of  molecules  or  exotic 
atomic  species  [8,9],  and  quantum  state  engineering  of 
spins  in  semiconductors  [10,11]  or  ro vibrational  states  in 
molecules  [12,13].  The  optical  frequency  comb  may  be¬ 
come  a  crucial  component  in  the  held  of  quantum  infor¬ 
mation  science,  where  complex  multilevel  quantum 
systems  must  be  controlled  with  great  precision  [14]. 

In  this  Letter,  we  report  the  use  of  an  optical  frequency 
comb  generated  from  an  ultrafast  mode-locked  laser  to 
efficiently  control  and  faithfully  entangle  two  trapped 
atomic  ion  qubits.  The  optical  pulse  train  drives  stimulated 
Raman  transitions  between  hyperhne  levels  [15,16],  ac¬ 
companied  by  qubit  state-dependent  momentum  kicks 
[17].  The  coherent  accumulation  of  these  pulses  generates 
particular  quantum  gate  operations  that  are  controlled 
through  the  phase  relationship  between  successive  pulses. 
This  precise  spectral  control  of  the  process  along  with  the 
large  optical  bandwidth  required  for  bridging  the  qubit 
frequency  splitting  forms  a  simple  method  for  controlling 
both  the  internal  electronic  and  external  motional  states  of 
trapped  ion  qubits,  and  may  be  extended  to  most  atomic 
species.  This  same  approach  can  be  applied  to  control 
larger  trapped  ion  crystals  with  more  advanced  pulse¬ 
shaping  techniques,  and  can  also  be  extended  to  a  strong 
pulse  regime  where  only  a  few  high-power  pulses  are 
needed  for  fast  quantum  gate  operations  in  trapped  ions 
[17-19]. 

High  fidelity  qubit  operations  through  Raman  transitions 
are  typically  achieved  by  phase-locking  frequency  compo¬ 
nents  separated  by  the  energy  difference  of  the  qubit  states. 
This  is  traditionally  accomplished  in  a  bottom-up  type  of 
approach  where  either  two  monochromatic  lasers  are  phase 
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locked  or  a  single  cw  laser  is  modulated  by  an  acousto¬ 
optic  (AO)  or  an  electro-optic  (EO)  modulator.  However, 
the  technical  demands  of  phase-locked  lasers  and  the  lim¬ 
ited  bandwidths  of  the  modulators  hinder  their  application 
to  experiments.  Here  we  exploit  the  large  bandwidth  of 
ultrafast  laser  pulses  in  a  simple  top-down  approach  toward 
bridging  large  frequency  gaps  and  controlling  complex 
atomic  systems.  By  starting  with  the  broad  bandwidth  of 
an  ultrafast  laser  pulse,  a  spectral  landscape  can  be 
sculpted  by  interference  from  sequential  pulses,  pulse 
shaping  and  frequency  shifting.  In  this  Letter,  we  start 
with  a  picosecond  pulse  and,  through  the  application  of 
many  pulses,  generate  a  frequency  comb  that  drives  Raman 
transitions  by  stimulating  absorption  from  one  comb  tooth 
and  stimulating  emission  into  another  comb  tooth  as  de¬ 
picted  in  Fig.  1.  Because  this  process  only  relies  on  the 
frequency  difference  between  comb  teeth,  their  absolute 
position  is  irrelevant  and  the  carrier-envelope  phase  does 
not  need  to  be  locked.  As  an  example  of  how  this  new 
technique  promises  to  ease  experimental  complexities, 
metastable-state  qubits  separated  in  frequency  by  a  tera¬ 
hertz  have  been  controlled  using  cw  lasers  phase-locked 
through  a  frequency  comb  [20],  but  might  be  controlled 
directly  with  a  100  fs  Ti: sapphire  pulsed  laser. 

At  a  fixed  point  in  space,  an  idealized  train  of  laser 
pulses  has  a  time-dependent  electric  field  that  can  be 
written  as 

N 

E{t)  =Y^f{t-  nT)e‘“‘‘,  (1) 

where  f{t)  is  the  pulse  envelope,  T  is  the  time  between 
successive  pulses  (repetition  rate  Vr  =  IjT),  N  is  the 
number  of  pulses  in  the  train  and  is  the  carrier  fre¬ 
quency  of  the  pulse.  For  simplicity,  any  pulse-to-pulse 
optical  phase  shift  is  ignored  since  the  offset  frequency 
in  the  comb  is  unimportant.  The  Fourier  transform  of 
Eq.  (1)  defines  a  frequency  comb  characterized  by  an 
envelope  f{(x))  =  S[/(0]  centered  around  the  optical  fre¬ 
quency  (x)c  and  teeth  separated  by  whose  individual 
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FIG.  1  (color  online).  The  Stokes  Raman  process  driven  by 
frequency  combs  is  shown  here  schematically.  An  atom  starting 
in  the  |i)  state  can  be  excited  to  a  virtual  level  by  absorbing  a 
photon  from  the  blue  comb  and  then  driven  to  the  |t)  state  by 
emitting  a  photon  into  the  red  comb.  Although  drawn  here  as  two 
different  combs,  if  the  pulsed  laser’s  repetition  rate  or  one  of  its 
harmonics  is  in  resonance  with  the  hyperfine  frequency,  the  ab¬ 
sorption  and  emission  can  both  be  stimulated  by  the  same  fre¬ 
quency  comb.  Because  of  the  even  spacing  of  the  frequency 
comb,  all  of  the  comb  teeth  contribute  through  different  virtual 
states  which  result  in  indistinguishable  paths  and  add  construc¬ 
tively. 


widths  scale  like  ^v^/N.  The  Raman  resonance  condition 
will  be  satisfied  when  a  harmonic  of  the  repetition  rate  is 
equal  to  the  hyperfine  qubit  splitting  mq,  implying  that  the 
parameter  q  =  mqIIttvji,  is  an  integer.  To  demonstrate 
coherent  control  with  a  pulse  train,  ^^^Yb^  ions  con¬ 
fined  in  a  linear  Paul  trap  are  used  to  encode  qubits  in 
the  ^5*^/2  hyperfine  clock  states  \F  =  0,  m/7  =  0)  =  |i)  and 
\F  =  1,  m/7  =  0)  =  It),  having  hyperfine  splitting 
(OqIItt  =  12.6428  GHz.  Eor  state  preparation  and  detec¬ 
tion  we  use  standard  Doppler  cooling,  optical  pumping, 
and  state-dependent  fluorescence  methods  on  the  811  THz 
^‘^1/2  ^^1/2  electronic  transition  [21].  The  frequency 

comb  is  produced  by  a  frequency  doubled  mode-locked 
Ti:sapphire  laser  at  a  carrier  frequency  of  802  THz,  de¬ 
tuned  by  A/27r  =  9  THz  from  the  electronic  transition. 
The  repetition  rate  of  the  laser  is  =  80.78  MHz,  with 
each  pulse  having  a  duration  of  r  ~  1  psec.  The  repetition 
rate  is  phase-locked  to  a  stable  microwave  oscillator  as 
shown  in  Pig.  2,  providing  a  ratio  of  hyperfine  splitting  to 
comb  spacing  of  ^  =  156.5.  An  EO  pulse  picker  is  used  to 
allow  the  passage  of  one  out  of  every  n  pulses,  decreasing 
the  comb  spacing  by  a  factor  of  n  and  permitting  integral 
values  of  q.  As  shown  in  Pig.  3,  when  n  =  2  (q  =  313  and 
=  40.39  MHz),  application  of  the  pulse  train  drives 
oscillations  between  the  qubit  states  of  a  single  ion. 
However,  when  n  =  3  (^  =  469.5  and  =  26.93  MHz), 
the  qubit  does  not  evolve. 

The  Rabi  frequency  of  these  oscillations  can  be  esti¬ 
mated  by  considering  the  Hamiltonian  resulting  from  an 


FIG.  2  (color  online).  Schematic  of  the  experimental  setup 
showing  the  paths  of  the  pulse  trains  emitted  by  a  mode-locked 
Ti:  sapphire  (Ti:sapph)  laser,  where  the  optical  pulses  are  fre¬ 
quency  shifted  by  AOs.  Single  qubit  rotations  only  require  a 
single-pulse  train,  but  to  address  the  motional  modes  the  pulse 
train  is  split  into  two  and  sent  through  AOs  to  tune  the  relative 
offset  of  the  two  combs.  We  lock  the  repetition  rate  (^^)  by  first 
detecting  vr  with  a  photodetector  (PD).  The  output  of  the  PD  is 
an  RF  frequency  comb  spaced  by  Vr.  We  bandpass  filter  (BP)  the 
RF  comb  at  12.685  GHz  and  then  mix  the  signal  with  a  local 
oscillator  (LO).  The  output  of  the  mixer  is  sent  into  a  feedback 
loop  (PID)  which  stabilizes  by  means  of  a  piezo  mounted  on 
one  of  the  laser  cavity  mirrors.  When  locked,  is  stable  to 
within  1  Hz  for  more  than  an  hour.  As  an  alternative,  instead  of 
locking  the  repetition  rate  of  the  pulsed  laser,  an  error  signal 
could  be  sent  to  one  of  the  AOs  to  use  the  relative  offset  of  the 
two  combs  to  compensate  for  a  change  in  the  comb  spacing. 


infinite  train  of  pulses.  After  adiabatically  eliminating  the 
excited  ^^1/2  state  and  performing  the  rotating-wave  ap¬ 
proximation,  the  resonant  Rabi  frequency  of  Raman  tran¬ 
sitions  between  the  qubit  states  is  given  by  a  sum  over  all 
spectral  components  of  the  comb  teeth  as  indicated  in 
Fig.  1  {h  =  1): 


l-q 


fin 


(x)qT 


)■ 


(2) 


where  fjL  is  the  dipole  matrix  element  between  the  ground 
and  excited  electronic  states,  and  q  is 

an  integer.  In  the  approximate  expression  above,  the  sum  is 
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FIG.  3  (color  online).  After  Doppler  cooling  and  optical  pump¬ 
ing  to  the  I  i)  state,  a  single-pulse  train  is  directed  onto  the  ion. 
When  the  ratio  of  qubit  splitting  to  pulse  repetition  rate,  is  an 
integer,  pairs  of  comb  teeth  can  drive  Raman  transitions  as 
shown  by  the  blue  circular  data  points.  However,  if  the  q 
parameter  is  a  half-integer,  the  qubit  remains  in  the  initial  state 
as  shown  by  the  red  square  data  points. 
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replaced  by  an  integral  and  each  pulse  is  described  by 
/(^)  =  V^r72£'osech(7rf/T)  with  r  T,  where  Oq  = 
(^'/^T)|yLt£’oP/A  =  sy^l2A  is  the  time-averaged  resonant 
Rabi  frequency  of  the  pulse  train,  yllir  =  20  MHz  is  the 
radiative  linewidth  of  the  ^  ^  is  the 

average  intensity  7=  Pj^c€qI2  J dt\f{t)\^  scaled  to  the 
^‘^1/2  ^^1/2  saturation  intensity.  Note  the  net  transition 

rate  is  suppressed  unless  the  single-pulse  bandwidth  is 
large  compared  to  the  hyperfine  frequency  (coqT  1),  in 
which  case  ft  ~  fto-  In  our  experiments,  coqT  ~  0.08.  Eor 
/sat  =  0.15  W/cm^,  the  data  shown  in  Eig.  3  is  consistent 
with  an  average  intensity  7  —  500  W/cm^. 

In  order  to  entangle  multiple  ions,  we  first  address  the 
motion  of  the  ion  by  resolving  motional  sideband  transi¬ 
tions.  As  depicted  in  Eig.  2,  two  perpendicular  pulse  trains 
with  wave  vector  difference  k  along  the  v  direction  of 
motion  are  directed  onto  the  ion.  Their  polarizations  are 
mutually  orthogonal  to  each  other  and  to  a  weak  magnetic 
field  that  defines  the  quantization  axis  [22].  To  ensure 
temporal  overlap,  the  path  lengths  of  the  two  beams  are 
matched  to  better  than  c/r  — 0.3  mm.  We  control  the 
spectral  beat  notes  between  the  combs  by  sending  both 
beams  through  AO  modulators  (driven  at  frequencies  Vi 
and  V2),  imparting  a  net  offset  frequency  of  AcoI27t  = 
V\  —  Vi  between  the  combs.  Eor  instance,  in  order  to  drive 
the  first  upper  or  lower  sideband  transition  we  set 
+  Aa)|  =  coq±  cOf,  with  j  an  integer  and  the 
trap  frequency.  In  order  to  see  how  the  sidebands  are 
spectrally  resolved,  we  consider  the  following 
Hamiltonian  of  a  single  ion  and  single  mode  of  harmonic 
motion  interacting  with  the  Raman  pulse  train: 

+  ^0 

T/eff  =  a  + 

^  n 

(3) 

where  =  ftT  is  the  change  in  the  Bloch  angle  due  to  a 
single  pulse,  is  the  Pauli-z  operator,  a±  are  raising  and 
lowering  operators,  x  is  the  v-position  operator  of  the 
trapped  ion,  and  a  are  the  raising  and  lowering  operators 
of  the  V  mode  of  harmonic  motion  and  the  q  parameter  has 
been  assumed  to  not  be  an  integer  or  half-integer.  In  the 
interaction  picture,  the  evolution  operator  after  N  pulses  is 
given  by  where 

_ -i 

— (4) 

and  Hq  =  a  +  \/2{(x)q  +  A(o)a^.  The  time  evolution 
operator  is  given  by, 

/  6  \ 

=  exp[- /HoA^r](  /  -  /  ^  X  )’  (5) 

^  «=0  ^ 

^  +  H.C.,  (6) 

where  D{a)  =  exp[afa^  —  Qf*a]  is  the  harmonic  oscillator 


displacement  operator  in  phase  space,  and  7]  =  k^lhl2m(o^ 
is  the  Lamb-Dicke  parameter.  In  the  Lamb-Dicke  regime, 

a)  +  1  1,  we  can  write  — 

1  +  turning  the  sum  in  Eq.  (5) 

into  a  geometric  series.  If,  for  example,  the  offset  fre¬ 
quency  between  the  combs  A  co  is  tuned  to  satisfy  the 
resonance  condition  for  the  red  sideband,  = 
{(Oq  +  a (o  —  (Of)T  =  277 j,  where  j  is  an  integer,  then  the 
sum  in  Eq.  (5)  is  approximately  given  by 


v-i 


I 


Qn  ^  iV 


smNd-,/2 

sin#^/2 


(V) 


The  coefficient  in  Eq.  (7)  is  the  same  as  the  field  amplitude 
created  by  a  diffraction  grating  of  N  slits,  whose  narrow 
peaks  have  an  amplitude  equal  to  N.  In  the  limit  co^T  I, 
the  other  terms  in  Eq.  (5)  that  drive  the  carrier  and  other 
sideband  transitions  can  be  neglected  when  N » 
(cOfT This  is  analogous  to  the  destructive  interference 
of  amplitudes  away  from  the  bright  peaks  in  a  diffraction 
grating.  Eor  cl>J27t  =  1.64  MHz,  T  =  12.4  ns,  and  rj  = 
0.1,  the  sidebands  are  well  resolved  when  N  »  80. 


AOl  Frequency  (MHz) 


Detuning  from  carrier  transition  (MHz) 


FIG.  4  (color  online),  (a)  Using  a  Raman  probe  duration  of 
80  yus,  (A  —  6500),  a  frequency  scan  of  AOl  shows  the  resolved 
carrier  and  motional  sideband  transitions  of  a  single  trapped  ion. 
The  transitions  are  labeled,  (An^,  Ariy),  to  indicate  the  change  in 
the  number  of  phonons  in  the  two  transverse  modes  that  accom¬ 
pany  a  spin  flip.  The  x  and  y  mode  splitting  is  controlled  by 
applying  biasing  voltages  to  the  trap  electrodes.  Unlabeled  peaks 
show  higher  order  sideband  transitions  and  transitions  to  other 
Zeeman  levels  due  to  imperfect  polarization  of  the  Raman 
beams,  (b)  Ground  state  cooling  of  the  motional  modes  via  a 
train  of  phase-coherent  ultrafast  pulses.  The  red  open-circle  data 
points  show  that  after  Doppler  cooling  and  optical  pumping,  both 
the  red  and  blue  sidebands  are  easily  driven.  The  blue  filled- 
circle  data  points  show  that  after  sideband  cooling,  the  ion  is 
close  to  the  motional  ground  state,  (n^y  <  0.03),  as  evidenced 
by  the  suppression  of  the  red-sideband  transition. 
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FIG.  5  (color  online).  The  parity  oscillation  that  is  used  to 
calculate  the  fidelity  of  the  spin  state  of  two  ions  with  respect  to 
the  maximally  entangled  state  |;^)  after  performing  the  entan¬ 
gling  gate.  The  phase  </>  of  the  analyzing  pulse  is  varied  by 
changing  the  relative  phase  of  the  rotation  pulses;  the  values  of  <p 
were  chosen  for  efficiency  in  determining  the  parity  contrast. 

Eor  many  applications  in  quantum  information,  the  mo¬ 
tional  modes  of  the  ion  must  be  cooled  and  initialized  to  a 
nearly  pure  state.  Eigure  4  shows  that  the  pulsed  laser  can 
also  be  used  to  carry  out  the  standard  techniques  of  side¬ 
band  cooling  [22]  to  prepare  the  ion  in  the  motional  ground 
state  with  near  unit  fidelity.  The  setup  also  easily  lends 
itself  to  implementing  a  two-qubit  entangling  gate  by 
applying  two  fields  whose  frequencies  are  symmetrically 
detuned  from  the  red  and  blue  sidebands  [23,24].  By 
simultaneously  applying  two  modulation  frequencies  to 
one  of  the  comb  AO  frequency  shifters,  we  create  two 
combs  in  one  of  the  beams.  When  these  combs  are  tuned 
to  drive  the  red  and  blue  sidebands  (in  conjunction  with  the 
third  frequency  comb  in  the  other  beam),  the  ion  experi¬ 
ences  a  spin-dependent  force  in  a  rotated  basis  as  described 
in  Ref.  [25].  Ideally,  when  the  fields  are  detuned  from  the 
sidebands  by  an  equal  and  opposite  amount  5  =  27/0,  a 
decoupling  of  the  motion  and  spin  occurs  at  gate  time  = 
Itt 1 8,  and  the  spin  state  evolve  to  the  maximally  entangled 
state  Ia)  =  lU)  +  In  the  experiment,  =  108  yas 

(V  ~  8700  pulses). 

The  fidelity  of  the  state  p  with  respect  to  \x)  is  found  by 
measuring  the  populations  pp  n  and  and  varying  the 
measurement  phase  angle  cf)  in  the  parity  signal  n((/>)  = 
T r[^(T^P pR^  where  R{4>)  is  a  7r/2  rota¬ 

tion  on  the  Bloch  sphere  with  phase  </>.  The  contrast  of  the 
parity  signal,  H^;,  is  used  to  calculate  the  fidelity  = 
(Ptftt  Pu,u)/^  +  nc/4-  The  measured  populations  of  |U) 
and  It!)  together  with  the  data  shown  in  Eig.  5  yield  a 
fidelity  =  0.86  ±  0.03,  thereby  signaling  that  the  two 
ions  are  entangled  after  the  application  of  the  pulse  train 
[26].  The  measured  fidelity  is  limited  primarily  by  an 
estimated  —10%  error  in  determining  the  two-ion  spin 
state  through  histogram  fitting.  This  error  stems  largely 
from  uneven  illumination  of  the  ions  during  detection, 
leading  to  a  bias  in  the  parity  signal  as  seen  in  Eig.  5. 
Errors  from  power  and  beam-steering  fluctuations  are  ex¬ 
pected  to  contribute  —4%  to  the  infidelity.  The  rate  of  off- 
resonant  spontaneous  photon  scattering  is  expected  to  be 


Rs  =  2yCl/A,  assuming  the  laser  pulse  envelope  f(t)  is 
smooth  and  has  transform-limited  bandwidth  —  l/r  A. 
The  probability  of  spontaneous  emission  during  the  entan¬ 
gling  gate  operation  is  thus  expected  to  be  R^tg  = 
liry  I  r/A  —  10“"^. 

We  have  demonstrated  full  control  and  entanglement  of 
two  atomic  qubits  using  an  optical  frequency  comb.  This 
work  represents  a  significant  simplification  over  current 
methods  for  optical  control  of  trapped  ion  qubits,  and  also 
points  the  way  toward  future  advances  with  higher  power 
laser  pulses.  Eor  example,  such  pulses  allow  much  larger 
detunings  from  resonance  and  a  suppression  of  decoher¬ 
ence  from  spontaneous  emission  [27]  while  maintaining 
gate  speed.  When  only  a  few  high-power  pulses  are  used  in 
ways  similar  to  the  experiment,  it  also  becomes  possible  to 
suppress  sources  of  motional  decoherence  through  the  use 
of  fast  entanglement  schemes  [18]. 
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